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Abstract 

The parametric equations of the plane curves determining the equihbrium 
shapes that a uniform inextensible elastic ring or tube could take subject to 
a uniform hydrostatic pressure are presented in an explicit analytic form. The 
determination of the equilibrium shape of such a structure corresponding to a 
given pressure is reduced to the solution of two transcendental equations. The 
shapes with points of contact and the corresponding (contact) pressures are 
determined by the solutions of three transcendental equations. The analytic 
results presented here confirm many of the previous numerical results on this 
subject but the results concerning the shapes with lines of contact reported 
up to now are revised. 
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1. Introduction 



In the present paper, the problem for determination of the equilibrium 
shapes of a circular inextensible elastic ring subject to a uniformly distributed 
external force that acts normally to the ring in the ring plane is addressed. 
This problem is also referred to as the stability problem or buckling of the 
circular shape of the ring and the other equilibrium shapes are called buckled 
It is also known (see, e.g., 0, O, Q]), that if a cylindrical elastic 
shell of circular cross section (i.e., a tube) is subject to a uniform external 
pressure, which is normal to its middle surface, then the typical cross section 
of the deformed tube takes the same shapes as the axis of a deformed elastic 
ring does provided that the latter is a simple curve (i.e., a curve without 
intersections). Therefore, here the term "ring" will be used to indicate both 
a proper ring and a tube. It should be noted also that in the majority of the 
works in this field, the distributed force acting on a ring is called pressure as 
in the case of a shell. Following this tradition, we will use the same term in 
the present study remembering that pressure means force per unit length in 
the case of a ring and force per unit area in the case of a shell. 

Maurice Levy [t] was the first who stated and studied the problem under 
consideration and reduced the determination of the foregoing equilibrium 
shapes in polar coordinates to two elliptic integrals for the arclength and 
polar angle regarded as functions of the squared radial coordinate. He found 
also several remarkable properties of the equilibrium ring shapes and obtained 
that if the pressure p is such that p < (9/4) (D/p^), where D and p are the 
ring bending rigidity and radius of the undeformed shape, respectively, then 
the ring possesses only the circular equilibrium shape. 
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Later on, Halphen [8| and Greenhill {q] derived exact solutions to this 
problem in terms of Weierstrass elliptic functions on the ground of com- 
plicated analyses of the properties of the aforementioned elliptic integrals. 
Halphen (see sl, p. 235]) found out that non-circular shapes with n > 2 axes 
of symmetry are possible only for pressures greater than p„ = {n'^ — 1){D/ p^). 
Halphen jsj] and Greenhill 9| presented also several examples of non-circular 
equihbrium ringshapes. It should be noted, however, that the exact solutions 



reported 



im rmgsi 
in 0,3, 



representing the polar angle as a function of the radius. 



appeared to^ intractable and many researchers continued searching exact 

151], while others used various approximations 




16| on the way to determine the equilibrium shapes of the ring. 



Carrier [ij] was the first who reconsidered the foregoing problem for the 
buckling of an elastic ring about half a century after the works by Levy, 
Halphen and Greenhill. He expressed the curvature of the deformed ring in 
terms of Jacobi cosine function 17] involving several unknown parameters to 
be determined by a system of algebraic equations. However, he succeeded to 
find approximate solutions to this system only for small deflections from the 
undeformed circular ring shape (see the exhaustive analysis provided recently 
11 1 who has criticized and developed Carrier's work jl|). 



by Adams 



Tadjbakhsh and Odeh [2[ studied the boundary-value problem describing 
the buckled shapes of the ring and the associated variational problem. They 
proved the existence of solutions to the boundary- value problem in the case of 
small deflections from the undeformed state and the existence of solutions of 
the associated variational problem (week solutions to the foregoing boundary- 
value problem) describing buckled shapes of an arbitrary deflection. 
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Watanabe and Takagi 12| thoroughly analyzed the variational problem 
for determination of the ring shapes stated by Tadjbakhsh and Odeh 
and obtained analytic expressions for the curvature of the ring (in terms of 
Jacobi elliptic functions and Carrier's parameters pj) and formulas for the 
slope angle at the points where the curvature has extrema. They also proved 
that non-circular shapes with n > 2 axes of symmetry exist for pressures 
greater than p„ = (n^ — l){D/p^), thus extending Halphen's result, and 
found out, moreover, that each such shape in unique. This paper completes 
the branch of analysis of the considered problem that is based more or less 
on the approach suggested more than half a century ago by Carrier jlj. 



In the recent papers 
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14 



15| . the present authors have studied the 



differential equation for the curvature of the ring with the aim to achieve 
an analytic description of the cylindrical equilibrium shapes of lipid bilayer 
membranes. The determination of the analytic solutions of this equation 
reported in 15| does not follow Carrier's approach jlj|. Instead of this, the 
explicit formulas for the curvature of the buckled shapes are obtained in 
forms similar to those suggested by Zhang 10| for lipid bilayer membranes 
and by Fukumoto 18| in the context of fluid mechanics. In [15], the para- 



metric equations of the directrices of the considered cylindrical surfaces are 
expressed in an explicit analytic form, the necessary and sufficient conditions 
for such a surface to be closed are derived and several sufficient conditions 
for its directrix to be simple or self-intersecting are given. 

The equilibrium shapes of closed planar elastic loops subject to the con- 
straints of fixed length and enclosed area are studied also in the works by 



Arreaga et al. 



19| . Capovilla et al. 



20| and Guven [21 1. 
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Flaherty et al. ^\ presented a numerical determination of the equilibrium 
shapes of elastic rings or tubes subject to uniform external pressure. They 
suggested a scenario for the evolution of the equilibrium shape as the pressure 
increases, however, some stages of this scenario are not confirmed here. 

Recently, Wu et al. [l6| presented approximate solutions of the considered 



problem by means of a special Fourier expansion of the curvature. They 
developed a linearized algebraic system for the unknown coefficients of this 
Fourier expansion and thus obtained an approximation of the ring shape. 

The aim of the present paper is to provide an analytic description of the 
equilibrium states of an elastic ring or tube subject to a uniform hydrostatic 
pressure going as far as possible, to develop efficient computational proce- 
dures completing this analytic description and to reexamine some of the most 
important results obtained in this field. 

The layout of the paper is as follows. The statement of the problem is 
given in Section 2. A concise derivation of the most important Levy's results 
is given in Section 3 where the equation for the ring curvature and the para- 
metric equations of the equilibrium shapes are derived as well. Using these 
results, in Section 4, the symmetry of the equilibrium shapes is justified. All 
periodic solutions to the equation for the ring curvature and the expressions 
for the corresponding slope angles are presented in explicit analytic form in 
Section 5. Two systems of transcendental equations are derived in Sections 6 
and 7 allowing to determine the equilibrium shapes corresponding to a given 
pressure and to calculate the values of the pressure at which the equilibrium 
shapes possess points of contact. The results concerning the equilibrium 



shapes with lines of contact obtained in 



J] are reexamined in Section 8. 
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2. Differential equations for the equilibrium states 

Let us consider a ring made of a homogeneous isotropic linearly elastic 
material and assume that it is represented by its middle axis. Suppose now 
that the ring is subject to a uniform external pressure p acting along the 
normal vector to its stress-free configuration. The analysis carried out in the 
present work is based on the following three assumptions: (z) the ring axis 
is inextensible; (ii) the pressure p preserves its magnitude and always acts 
as an external uniformly distributed force along the inward normal vector 
to the deformed ring axis, i.e., it is a uniform (simple) hydrostatic pressure; 
{in) the deformed ring axis is a regular closed plane curve F parametrized 
by its arclength s. 

Next, let the curve T be given by means of the coordinates x{s), y{s) 
of its position vector r(s) with respect to a certain rectangular Cartesian 
coordinate frame in the Euclidean plane, i.e., r(s) = x{s) i + y{s)j, where i 
and j are the unit vectors along the coordinate and y, respectively. 

Consequently, the unit tangent vector t(s) and the unit inward normal vector 
n(s) to the curve F are given as follows 

t{s) = x'{s)i + y'{s)], n{s) = -y'{s)i + x'{s)]. (1) 

Here and throughout this paper primes denote derivatives with respect to 
the arclength s. Recall that the foregoing unit tangent and normal vectors 
are related to the curvature k(s) of the curve F through the Frenet-Serret 
formulas 

t'(s) = K{s)n{s), n'(s) = -K{s)t{s). (2) 
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Finally, let M{s), N{s) and Q{s) denote the bending moment and the 
components of the stress resultant force F(s) along the tangent and normal 
vectors to the curve F, i.e., 

F(s) = Ar(s)t(s) + g(s)n(s). (3) 

Then, the particular constitutive equation relating the moment M{s) with 
the curvature k{s) and the form of the stress-free configuration of the ring 
together with the system of differential equations 

F\s) = -pn{s), (4) 

M'{s) = -F{s)-n{s), (5) 

representing the local balances of the force and moment, respectively, in 
accordance with the assumption (m), and the closure conditions following 
from the assumption {Hi), which, without loss of generality, may be written 
in the form 

r(L) = r(0), t(L)=t(0), (6) 

where L is the length of the deformed ring, determine entirely the equilib- 
rium state of the ring under consideration (see, e.g., jsl, Ch. 4)]). Here and 
throughout this paper the dot stands for dot (scalar) product of two vectors. 
Let us remark also that using Eq. ([3]) and the Frenet-Serret formulas (|2]) one 
can represent the system of differential equations (jl]), ([5]) in the scalar form 

N'{s) = Qis)K{s), (7) 
Q'is) = -N{s)K{s)-p, (8) 
M'{s) = -Q{s). (9) 



3. Parametric equations for the equilibrium shapes 

The aim of this Section is to prove the most important facts concerning 
the problem under consideration established by Maurice Levy in his memoir 
[3] and to derive the parametric equations of the equilibrium ring shapes. 

Using the expression for the normal vector, see Eqs. ([T]), one can integrate 
the equilibrium condition and express the force vector in the form 

F{s) = py{s)i- px{s)}, (10) 

omitting the constant of integration since it is always possible to eliminate it 
by choosing the origin of the coordinate frame at a certain privileged point 
which, actually, is the one that Levy called "centre of the elastic forces", 
cf. 0, 1° (p- 9)]. Eq. (ITOl) implies that at each point of the deformed ring 
configuration the magnitude of the force vector F{s) = |F(s)| is proportional 
to the magnitude of the position vector r{s) = |r(s)|, the magnitude of the 
pressure \p\ being the coefficient of proportionality (cf. 3, 2° (p. 9)]), i.e., 

F{s) = \p\r{s). (11) 

Next, taking the dot product of both sides of Eq. (fTO!) with the normal vector 
one gets, bearing in mind the second one of equations ([1]), the relation 

F(s) ■ n(s) = -p[x{s)x'{s) + y{s)y'{s)], 

which allows to integrate the balance of moment equation ([5]) and to obtain 

M{s) = ^{r\s) + c), (12) 

n 

where c is an arbitrary constant of integration, cf. [7, 3° (p. 9)]. It is 
noteworthy that the relations ffTOj) - f|T2|) hold regardless of the particular 
material properties of the ring and the form of its stress-free configuration. 
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In terms of the slope angle !f{s) one has the expressions 

x'{s) = cos ip{s), y'{s) = sirnp{s), (13) 
<s) = ^'(s), (14) 

and using Eqs. f|T3l) can rewrite Eqs. ([1]) in the form 

t(s) = cosv9(s) i + sin(y9(s) j, (15) 

n(s) = — sin(y9(s) i + cosv9(s) j. (16) 

Then, combining Eqs. ([3]), ( fTOl) . ( ITSi) and ( fT6l) one obtains the parametric 
equations of the deformed ring shape in the form 

x(s) = — Q(s) cos(p(s) N(s) sirnpis), 

p p 

y{s) = — Qis) sm(f{s) H — N{s) cosip{s). (17) 
p p 



Evidently, in view of Eq. ( ITSll . the second one of the closure conditions (EI) 
implies that the rotation number of the deformed ring axis is 2m7r, where m 
is an integer, i.e., 

<^{L) = <^{0) + 2m7r, (18) 

whereas the first one of them transforms, on account of Eqs. (fT2ll . (fT7I) and 
(fT8|) . into the obvious conditions for periodicity of the forces and moment 

Nil) = N{0), Q{L) = g(0), M(L) = M(0). (19) 



Again, neither the form of the parametric equations flTTl) nor the forms of 
the boundary conditions f|T8l) and f|T9l) depend on the particular material 
properties or the stress-free configuration of the ring. 
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Let us now assume that the constitutive equation of the ring is 

M{s) = D{k{s) - K°), (20) 

where D is its bending rigidity and k° = 1/p is the curvature of its stress-free 
configuration, which is supposed to be a circle of radius p. Using Eq. fl2U]) 
and the Frenet-Serret formulas ([2]) we can immediately integrate the system 
of differential equations ((7]) - (|9]) obtaining the following expressions for the 
tangent and normal components of the stress resultant force F(s) 

N{s) = -^{k\s) - 2/i), Q{s) = -Dk'{s), (21) 

and the single ordinary differential equation for the ring curvature 

n"{s) + ^K^is) - pk{s) - a = 0, (22) 

where a = p/D and /i is an arbitrary constant of integration. On the other 
hand, combining equations (!20!) . (ITT]) and (fT2l) we obtain the relation 

N\s) + Q\s) = 2pD{K{s) - K°) - p^c, 
which, in view of Eqs. f l2T|) . implies 

K\sr = PiK{s)), (23) 
where P{n) is a fourth-order polynomial of the curvature k of the form 

P{k) = ~^K^ + fXK^ + 2aK + e (24) 

whose free term e = —2o'k° — a'^c — p^ incorporates all the constants of 
integrations introduced so far. Actually, Eq. fj23l) is a first integral of Eq. fj22|) 
(see |15, Sec. 2] for more details). In this context, e is viewed as an arbitrary 
constant of integration. 
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Each sufficiently smooth real-valued solution k,{s) of an equation of form 
fl23|) corresponding to a certain triple of given values of the parameters /i, e 
and 0" 7^ generates, up to a rigid motion in the plane, a unique plane curve 
r of curvature k{s). The components of the position vector of this curve can 
by expressed in the form 

x{s) = —k'{s) cosip{s) + — — 2/i) sin(y9(s), 
a 2a 

y{s) = —k'{s) sm(f{s) — —{^'^{s) — 2/i) cosip{s), (25) 
a 2a 

obtained by substituting Eqs. fl2T|) in the general formulas f|T7|) . However, 
the parametric equations fl25l) describe a shape that a ring of bending rigidity 
D could take being subject to pressure p = aD if and only if the regarded 
solution K,{s) of the respective equation of form fl2^ is such that the closure 
conditions (fT8|) and (fT9|) hold for L = 27rp; note that the latter equality 
follows from the assumption (i). If this is the case, then the respective 
solution k{s), its first derivative k'(s) and its indefinite integral 

^p(s) = j K(s)ds, (26) 

cf. Eq. f|T^ . determine entirely the equilibrium state of the pressurized ring. 
Indeed, the shape of the ring is determined explicitly by Eqs. ( I25l) and the 
values of the moment and forces acting along the ring axes are given by the 
formulas ( 120|) and ( 12T|) . respectively. 

It should be remarked that each such solution k{s) is necessarily a periodic 
function with period L = 27rp, due to the condition fll9p . and if T is its least 
period, then L = nT, where n is a positive integer. Since ^p{nT) = ruplT), 
as follows by formula fl26|) . the closure condition f|T8|) takes the form 

^{T) = -y, + 27 

n n 
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4. Symmetry of the ring shapes 

The symmetry of the ring shapes is discussed, usually without going into 
much detail, by almost all authors that contributed to the solution of the 
considered problem. A concise justification of this property is given below. 

Let r be a curve whose curvature n{s) is a periodic solution of Eq. f l23|) 
with least period T. For each integer i denote by T^^^ and T^'^^ the parts of 
the curve F corresponding to s G [Tj — t, Tj] and s G [Tj, Tj + 1], respectively, 
where Ti = i{T/2) and t G M. Since the relations 

t(f,) ■ r(f, + s) = -t(f,) ■ r(f, - s) (28) 
n(f,) ■ r(f, + s) = n(r,) ■ r(f, - s) (29) 

hold for each couple of numbers i G N and s G M, the curves F^^^ ^ and F^-^^ 
are symmetric with respect to the axis directed by the normal vectors n(Tj). 
Fig. 1 illustrates the simplest case i = 1, t = Ti. 



y 




X 



Figure 1: Position vectors r{T/2 ± s), tangent t{T/2) and normal n{T/2) vectors to a 
curve; the dashed Hue represents the axis of symmetry. 
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It should be noted that relations and ([21]) are easily verified using 
the parametric equations (|25|1 and the formulas 

K{fi + s) = K{fi-s), K'{fi + S) = -K'{f,-S) (30) 

and 

ififi + s)= ififi) + ifis) (31) 

that follow from k(— s) = k(s) and the definition (l26l) of the slope angle f{s). 

In other words, the curve F can be thought of as generated by successive 
refiections of its part t[^ \ t = Ti, about the axes directed by the normal 
vectors n(Tj), i G N. Apparently, if the curve tI^^^ UF^^'' get closed smoothly 
for some positive integer n, which may happened only for t = Tn, then the 
closed curve F = F^^"* U '^nt^ has n distinct axes of symmetry and is said to 
have n-fold symmetry or to be of n mode. 



5. Analytic expressions for the ring curvatures and slope angles 



As it was mentioned in the Introduction, Eq. fl23|) was studied by the 
present authors with the aim to describe the cylindrical equilibrium shapes of 
fiuid lipid bilayer membranes (see 
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15|). Explicit analytic expressions 



for all periodic solutions of Eq. (!23|) and for the corresponding slope angles 
( 126|) are presented in 15|. These results are outlined in this Section. 

Depending on the values of the parameters a, fi and e, there exist two 
cases in which the polynomial P{n) attains positive values and hence Eq. f l23p 
has real-valued solutions: (I) the polynomial -P(k) has two simple real roots 
a, (3 E M., a < (3, and a pair of complex conjugate roots 7,5 G C, 5 = 7; 
(II) the polynomial -P(/t) has four simple real roots a</3<7<5G]R. In 
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the first case, the polynomial -P(k) is nonnegative in the interval a < k, < P, 
while in the second one, it is nonnegative in the intervals a < k < /3 and 
7 < K < 5. 

It should be noted that the roots a,/3,7 and 6 of the polynomial P{n) 
can be expressed explicitly through its coefficients /i, a and e and vice versa. 
Indeed, after some standard algebraic manipulations (see, e.g., |23|), one can 
find the following expressions for the roots of the polynomial P{k) 



-y/u ± J 2fl - 2(7^^ - 00, y/uj ± J 2fi + 2(7^^ 



where 



(2/x + C)^-2^3£ 
6C 



C = { 2^3 (3V2 + ^)- 23/. + 3%) 



X = 2'3£ 



2 o2 2 

— 6 fia 



3(T^ (2^/i^ - 3^2) . 



Then, one can denote properly each of the above expressions for the roots in 
accordance with the notation introduced in the cases (I) and (II), respectively. 
Simultaneously, by Vieta's formulas one obtains 



a + /3 + 7 + (5 = 0. 



(32) 



due to the absence of a term with in the polynomial -P(k), see Eq. (12^ . 
and 



/i = ^ (a^ + /J^ + 7^ + a/3 + a7 + /37) 
(T = -^{a + f3) (a + 7) (/3 + 7) , 
e = ^a(3-f{a + (3 +7). 



(33) 
(34) 
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5.1. Case (I) 

Let the parameters /i, a and e be such that the polynomial P{k) has 
roots as in case (I), namely, two of them are real (a < /3) and the other two 
constitute a complex conjugate pair which, in view of relation f l32p . can be 
written in the form 

a + 8. a + 13 . 
7 = Y~ 2 ^ ' 

where is a nonnegative real number. In this case, Eq. f |23p has periodic 

solutions if ?7 7^ or ?7 = and (3a; + /3)(q; + 3/3) > 0, see [l5|. Theorem 1]. 



Let rj ^ and hence the roots of the polynomial P{k,) are simple. Denote 



...ivii. ,,.^i-'-i±ss^±mi±m, (36) 



where 

A = ^V + (3« + /3)', B = ^J4r]^ + {a + 3f3f. (37) 

Evidently, A > 0, B > 0, Xi > and < A;i < L In this case, each solution 

of Eq. f l23|) can be expressed by the function 

_ {A/3 + Ba) - {AI3 - Ba) cn(Aig, feQ 

{A + B)-{A-B)cn{\,sM) ' ^ ^ 

which takes real values for each s G M and is periodic with least period 

Ti = (4/Ai)K(A;i) due to the periodicity of the Jacobi function cn(AiS, /ci). 
Here, K(-) denotes the complete elliptic integral of the first kind. The corre- 
sponding slope angle ( 126|) can be written in the form 

AI3 - Ba 



^1 (s) 



A-B 

+ arctan^/^.^"^^^^'^^^ 



2\i^kl + C dn(AiS,A;i 
+ ^-^^^^^^^ni-C,^mi\,s,h),h), (39) 

15 



where n(-, - , ■) denotes the incomplete elhptic integral of the third kind and 
C=<^. (40) 

4:AB ^ ' 

Now, let 7] = and (3a + [3){a + 3/3) > 0. Then, the polynomial -P(k) 
has one double and two simple real roots. The curvature and the slope angle 
are expressed in terms of elementary functions as follows 

{A(3 + Ba) - {AI3 - Ba) cos(Ais) 



1^2 [S] 



{A + B)- (A-5)cos(Ais) 



, , AP-Ba 8(a-/3) [A /I 

[s] = s + arctan I y ;g tan I -Ais 

where Ai, A and B remain defined by formulas (!36|) and (1371) . respectively. 
5.2. Case (II) 

Let the parameters /x, a and e be such that the polynomial -P(fi;) has roots 
as in case (II), that isa</3<7<5G]R. Denote 



Since 7 — a>/3 — q;>0 and 5 — /3 > 5 — 7 > 0, it is seen that A2 > and 
< ^2 < 1. In this case, each solution of Eq. ( l23l) can be expressed by one 
of the following functions 

{6-a){5-(3) 



^4(5) = /3 + 



{5- (3) + i(3-a)sn\X2S,k2)' 

(7-/3)(^-/3) 
(5-/3)-(<5-7)sn2(A2S,fc2)' 



see 
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, Theorem 2]. The functions n^is) and ^4(3) take real values for each 



s G M and are periodic with least period T2 = (2/A2)K(/c2) because of the 
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periodicity of the function sn^(A2S, ^2)- Their indefinite integrals (12^ can be 
written as 

^93(5) = 5s - — ^, am(A2S, ^2), h 



(3 — ^ f 5 — 
ipi{s) = (5s r — n -, am(A2S, ^2), k2 , 



A2 \(3-6 

'-7 -pr f 

A2 

respectively. 

It should be remarked, that the indefinite integrals ^Pj{s), j = 1, ... ,4, 
of the foregoing solutions i^j{s) of Eq. fl23|) are chosen so that v^j(O) = 0. 
Moreover, /tj(0) always coincides with a certain root of the polynomial -P(k) 
(actually, Kj{0) = a for j = 1,2,3 and ^4(0) = 7) and hence Kj(0) = 0, 
according to Eq. (!23|) . It is important to note also that the functions K,j{s) 
are strictly increasing for s G [0, T/2] where T is the respective least period. 

6. Determination of the equilibrium shapes 

In the present study, our primary interest is in the determination of the 
equilibrium ring shapes that are curves without intersections, i.e., simple 
curves. It is established in jl5| that only the solutions of Eq. f l2^ that fall 
under the case (I) with r] ^ and a + /3 7^ may give rise to simple curves. 
Therefore, hereafter we will restrict our analysis to the regular closed curves 
r„ of curvatures k,{s) = k,i{s) given by formula (!38|) . which have n axis of 
symmetry and meet all the necessary conditions for that to be simple. 

Thus, the closure condition (1271) for such a curve r„ reads 



^1 m) = ±-, (41) 

n 
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as ipi{0) = (see the remark at the end of Section 5), n > 2 due to the four 



vertex theorem (see, e.g., 2J]) and m = ±1 since the rotation number of a 
simple regular closed curve must be ±2tt, see j25|. Note, however, that there 
exist regular closed curve with rotation number ±2tt, which are not simple. 

Then, substituting the expression Ti = {4/ Xi)K{ki) for the least period of 
the solutions (138|) of Eq. ( 123|) in the general formula (!39|) for the corresponding 
slope angle ipi{s) one can rewrite the closure condition (14T]) in the form 

^-^±^^^^Ui-C, k,) + 4^Kik,) = ±^ . (42) 
2Xi{A-B) ^ ' Xi{A-B) ^ 2n ^ ^ 

Finally, substituting the same expression for the period Ti in the relation 

L = nTi in order to take into account that the length of the ring L is fixed 

and does not change upon deformation, see assumption (z), one obtains 

1k(.,) = iL . (43) 

after setting for simplicity, without loss of generality, L = 2tt, i.e., k,° = p = 1. 

The left-hand sides of Eqs. and fl43p are functions of the parameters 
a, (3 and t], see Eqs. fl36l) . fl37|) and PUI) . However, for the aims of the present 
study it is convenient to express the parameters a and f3, using formulas (IMj) 
and fl5^ . in terms of the positive parameters a, 1] and q as follows 

Thus, given an integer n > 2 and a pressure p by means of the parameter 
cr (called hereafter simply "pressure") the problem for the determination of 
the foregoing equilibrium shapes of the ring corresponding to this pressure 
is reduced to the computation of the solutions rj and q of the transcendental 
equations (|12D and fH5D . 
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It is important to notice that this problem has no nontrivial solution 



12 



if < o" < (Tbn and has a unique nontrivial solution if ct > o"b„, see 
Theorem 2]. Here, abn = — 1 is the so-called buckling pressure and by a 
trivial solution we mean the one, which corresponds to the ring shape that 
is a circle of radius p = 1. 

Here, the transcendental equations (H2|) and ( H3|) corresponding to given 
integer n > 2 and pressure a > abn are solved numerically in two steps 
using Mathematica®. First, the two curves in the (77,5) plane defined by 
Eqs. fH2|) and fH3|) are plotted using the routine ContourPlot in order to 
identify roughly the values of the coordinates of their intersection point. 
Then, these values are put as starting values in the routine FindRoot, which 
is employed to obtain the solutions r] and q of the system of equations fH2]) . 
( H3|) with a sufficient accuracy. Once such a solution is determined, formulas 
( H4|) . (!38|) . ( 139|) . fHOj) and the parametric equations ( 125|) allow to depict the 
corresponding equilibrium ring shapes using the routine ParametricPlot. 



Three examples of such shapes, which confirm the results presented in 
are given in Fig. 2. 





Figure 2: Equilibrium ring shapes corresponding to: (a) a = 4.75 (2-fold symmetry); (b) 
a = 16.25 (3-fold symmetry); (c) cr ~ 35.25 (4-fold symmetry). 
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7. Equilibrium shapes with points of contact 

In 2|, y], it is established that for each mode n = 2, 3, 4 there is a value 
of the pressure a, called contact pressure and denoted by den, at which some 
points of the respective buckled ring shape of n-fold symmetry come into 
contact. In the aforementioned works, it is also observed that if the applied 
pressure a is such that Ohn < o < Ocn-, then the corresponding buckled 
shape of n mode is simple. It should be noted, that the values for the 



contact pressures reported 



in 0,14 



^ are obtained solving numerically a rather 
complicated nonlinear boundary-value problem. 

In the present study, the determination of the non-circular equilibrium 
ring shapes with points of contact and the respective contact pressures is 
reexamined being reduced to the computation of the common solutions of 
the transcendental equations fH2l) . fH3l) and one more algebraic, when n = 2, 
or transcendental, when > 2, equation in the way to be described below. 

Proceeding to the examination of this problem, let us first clarify, slightly 
extending the definition used in ^, that an n-mode equilibrium ring shape 
r„ is said to have a point of contact if it is not self-intersecting, but there 
is at least one couple of values S\ and S2 of the archlength s such that 
< si < S2 < and 

r(s2) = r(si), t(s2) = ±t(si). (45) 

This means that at the point of contact r(s2) = the curve r„ is tangent 



to itself. Such a double point on a curve is called a cusp or tacnode, see [26|. 
The objective now is to reformulate the above conditions in a form suitable 
for the developing of an efficient procedure for computation of the contact 
pressures corresponding to the foregoing equilibrium ring shapes. 
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For that purpose, it is convenient to use the relations 

<s) = ^-As) - ^ (46) 

and 

r{s) ■ t{s) = -k'{s), r{s) ■ n{s) = -^{k\s) - 2/i), 

which follow from Eqs. (123|) - (125|) and allow, taking into account Eqs. (fT5i) 
and fll6p . conditions fl45p to be cast in the form 

Ki(s2) = Kl(Si), fi'l{s2) = ±fi[{Si), (47) 

Kl{s2)-2fi = ±{Kl{s,)-2fi), (48) 

and 

<fi{s2) =V>iisi) + 2h, (49) 

if the sign in the second one of Eqs. (HSI) is plus, or 

</'i(s2) =V^i(si) + (2/ + l)7r, (50) 

if the foregoing sign is minus. Here, / is an integer. 

In the case when t(s2) = t(si), conditions (H7|) imply S2 = Si +jTi, where 
j is a positive integer and hence, according to Eqs. fl26l) and fHTl) . we have 

2tt 

V9i(s2) = V5l(Si) ± j— ■ (51) 



Now, combining Eqs. ( 1491) and (jSTj) we obtain j = ±ln, which means that 
S2 = Si ± InTi = Si ± IL and therefore the assumptions t(s2) = t(si) and 
< Si < S2 < L turn out to be incompatible. That is to say that a non- 
circular equilibrium ring shape can not have cusps of this type. 
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In the remaining case t(s2) = — t(si), conditions ( H7|) and (HSj) read 

Ki{s2) = KliSi), k[{s2) = -k[{si), (52) 

kI{s2) -2fi = K,l{si) -2fi = 0. (53) 

Now, taking into account the symmetry of the considered type of ring shapes, 
the invariance of Eq. (!22|) under the translation of its independent variable s 
and the particular properties of its solutions reflected by the formulas 
without loss of generality one may assume that si G [0,Ti/2] and can easily 
enumerate all the values S2 of the arclength for which conditions fl52l) hold 
and the part of the curve r„ corresponding to s G [si, S2] may have double 
points of the considered type only, namely: S2 = Ti — si and S2 = nTi — Si. 
Hence, in view of Eqs. (126|) and ( 14T|) . condition (150|) takes the form 

^,{s,) = ±--^-ln, (54) 
n 2 

if S2 = Ti - si, or 

71" 

ipi{si) = ±n - - - In, (55) 
if S2 = nTi — Si. 

Next, suppose that the curve r„ is such that condition ( l53|l holds for some 
s G [0, Ti/2]. Actually, this means that in this interval there must be exactly 
two values o f t he arclength for which condition (1531) is fulfilled, otherwise, 
by virtue of p^. Theorem 3], the considered curve r„ is self- intersecting that 
contradicts to the accepted definition for curves with points of contact. Let 
us denote these values of the arclength by s]^ and s'l and assume that 



^i('5i ) = — V 2/i, Ki(sJ') = A/2/i . (56) 
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According to Eqs. there is some Sq G [0,Ti/2] for which k(so) = 

and therefore, in the hght of the above considerations, Si < sq < and 
'^i(O) < —V^Ji < < -\/2/i < Ki(Ti/2) because the curvature ni{s) is strictly 
increasing for s G [0,Ti/2] (see the remarks at the end of Section 5). Hence, 

kI{0) - 2/i > 0, /«?(Ti/2) - 2/i > . (57) 

Taking into account that k,'i{0) = n[{Ti/2) = and v'i(O) = 0, the first of 
parametric equations fl25|l imphes 

x(Ti/2) = ^ (k?(Ti/2) - 2/i) sin(^i(Ti/2). (58) 

The property fl3T]) of the slope angle allows the closure condition f HT]) to be 
recast in the form 

V9i(Ti/2) = (59) 

which, in view of the second of inequalities (1571) and the fact that a > 0, 
means that the sign of x(Ti/2) coincides with that of the right hand side 
of Eq. ( !59l) . If the latter sign is minus, then x(Ti/2) < as implied by 
Eq. ( l58l) and the second of inequalities ( 1571) . On the other hand, x(0) = 
and x'{s) = cos(f{s) meaning that x{s) attains positive values in the interval 
(0,Ti/2). Therefore, the respective curve r„ is self-intersecting in this case 
because it intersects one of its symmetry axis - the one directed by the normal 
vector n(0) = — j. In this way, we arrive at the conclusion that a contact 
point may occur only if the sign of the right hand side of closure condition 
Eq. ([59]) is plus. 
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It is clear that the slope angle v?i(s) has a minimum at s = Sq since 
'^[{sq) = /«i(so) = and iPi^so) = k[{sq) > 0. Actually, this is the only local 
extremum of this function in the interval (0,Ti/2) since Sq is the only point 
in the foregoing interval where the curvature is equal to zero. Therefore, 
(pi{s) < n/n for s G (0,Ti/2) since ^^(O) = and {p(Ti/2) = n/n. 

One can show also that ip{s) > —n for s G (0, Ti/2). Indeed, if ip{s) = — vr 
for some s G (0,Ti/2), then 

, _ x{s) cos (fjs) +y{s) siny?(g) _ _x(£) 
r[s) r[s) 

since x{s) > for each s in the considered interval. However, this inequality 
contradicts the fact that r(s) is strictly increasing for s G (0, Ti/2) as follows 
from Eq. (H6l) . and hence, ^p{s) > — tt in this interval. 

In view of the inequalities — vr < ip{s) < n/n, conditions (IMll and ( l55ll 
take the form 

77 71" 71" 

<^i(si) = ---, = --, 

respectively. Apparently, ^[{si) < and hence fi{s) is decreasing in the 
neighbourhood of s = while (y9']^(S]'") > and hence (pi{s) is increasing at 
s = s~[ . Therefore, contact points (if exist) should be such that 

IT TT , . 

Msl) = - - 2 ' (60) 

or 

^list) = -f , (61) 
otherwise the respective curve r„ is self-intersecting. 
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Finally, taking the inverse 



s = 



— F I arccos 



Af3 + Ba- k{A + B) 
A(3-Ba- n{A - B) 



,ki 



) 



of the function k, = which is readily achieved by Eq. ( 138|) and well 

defined for each s G [0,Ti/2], one, bearing in mind that s^,s^ G [0,Ti/2] 
and Eqs. obtains 



where F(-, ■) is the incomplete elliptic integral of the first kind. 

Thus, two triples of transcendental equations and or 

arise for the determination of the n-mode equilibrium ring shapes with points 
of contact. In both cases, the respective transcendental equations involve as 
unknowns only the four parameters a, t], q and n since the archlengths sj~ and 
sf, which may correspond to points of contact are determined explicitly in 
terms of these parameters by formulas fl62l) and (!63|) , and the same holds true 
for the left hand sides of Eqs. ( 160|) and ( 16T|) in view of the general expression 
( 13 9 p for the slope angle. Of course, one should remember that the sign of the 
right hand side of Eq. f H2|) is plus in this context. 

To summarize: given an integer n > 2, each solution of any one of the 
aforementioned two triples of transcendental equations gives the value of the 
contact pressure acn and the values of the parameters 1] and q determining in 
this way an equilibrium ring shape of n-fold symmetry with points of contact. 




(62) 




(63) 
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Solving numerically the foregoing two systems of transcendental equations 
using the routine FindRoot in Mathematica® we have found that the system 
consisting of Eqs. f l42|) . f l43|) and fl6T]) does not have solutions for 2 < n < 15, 
but the system of equations fj42l) . fH3|) and fl60|) has a unique solution for 
each such mode n. Our conjecture is that this happens for all modes. 

The obtained values for the respective contact pressures cxcn are presented 
in Table 1. The equilibrium ring shapes with points of contact corresponding 
to the contact pressures crc2, cfcz and are depicted in Fig. 3. 

Table 1. Contact pressures. 



n 


2 


3 


4 


5 


6 


7 


8 




5.247 


21.650 


51.844 


97.834 


161.077 


242.682 


343.517 


n 


9 


10 


11 


12 


13 


14 


15 




464.276 


605.522 


767.719 


951.253 


1156.450 


1383.580 


1632.890 



It is worth noting that in the special case n = 2 the transcendental 
equation ( l60i) may be replaced by the algebraic relation 

.^^Jl±fl. (64) 

16g ^ ' 

Indeed, in this case, expression fl60|) simplifies to ipi{s]^) = meaning that 
= 0. In fact, there exist two values of the archlength within the interval 
[0,Ti/2] in which (pi{s) = but the other one is necessarily greater than sq 
and therefore it is disregarded. Then, the first of Eqs. (156|) reads — 2/i = 0. 
Substituting here expressions fH^ i and ( 133|) for the root a and the parameter 
/i, respectively, and accounting for Eqs. and (H^ 9 one obtains Eq. 
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Figure 3: Ring shapes with points of contact: (a) a = 5.247; (b) o = 21.65; (c) a = 51.844. 



The values of contact jpressures cxcn obtained for n = 2, 3, 4 confirm exactly 
the results presented in [A, formulas (2.11)], see also Fig. 3, but the latter are 
interpreted therein as the lowest values of the pressures at which an isolated 
point of contact occurs. In , it is claimed that for each mode n beyond the 
contact pressure cXcn there exists a continuous range of pressures, from acn up 
to a certain pressure denoted by Uon, such that for each < cr < cron the 
respective ring shape exhibits contacts at isolated points only. The pressure 
(Ton is set in |4] to be the one for which the curvature at the corresponding 
contact point is zero. If so, however, then in view of Eqs. (66) fi = and 

= = sq. Therefore, the equation n(s) — 2u = has exactly one solution 



for s G [0,Ti/2] and hence, according to [l5|. Theorem 3], the corresponding 



ring shape is self-intersecting. Thus, in this respect the results of Flaherty 
et al. |4] turn out to be inaccurate in spite of the fact that they are widely 
accepted and even confirmed numerically by other authors (see, e.g., j^])- 

Actually, for all modes n in the range 2 < n < 15 our computations based 
on the procedure described in Section 6 show that if the applied pressure a 
is such that (Tf,„ < a < a cm then the corresponding buckled shape of n mode 
is simple, while for a > acn this shape always has points of self-intersection. 
Our conjecture is that this behaviour is inherent to all modes. 
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For n = 2, 3, 4, Flaherty et al. ^\ affirm that for each pressure a such that 
5.247 <a < 10.34, 21.65 <a < 81.81 or 51.84 <a < 207.2, respectively, the 
corresponding equilibrium ring shapes have only isolated points of contact 
without being self-intersecting. Our results presented in Figs. 4, 5 and 6 show 
that this is not the case. Let us recall that for each a > a hn t he corresponding 
equilibrium ring shape of n-fold symmetry is unique, see 12|, Theorem 2]. 





Figure 4: Ring shapes corresponding to: (a) a — 6.48; (b) a = 9.24; (c) a = 10.34. 






Figure 5: Ring shapes corresponding to: (a) a = 28.56; (b) a = 56.09; (c) a = 81.81. 






Figure 6: Ring shapes corresponding to: (a) a — 70.7; (b) a = 140; (c) a — 207.2. 
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8. Equilibrium shapes with hues (areas) of contact 

Apparently, for rings of finite thickness self-intersecting shapes are not 
possible because they are not planar, nevertheless for a very thin ring such 
a shape may be considered as a good approximation of its equilibrium state. 
A tube evidently can not take a self-intersecting shape, but there is a good 
reason to expect that tubes subject to sufficiently high pressure posses equi- 
librium shapes with areas of contact (lines of contact of their cross sections). 

Flaherty et al. [J] suggest similarity transformations to be used for the 
determination of such shapes, but realize this idea in a very complicated way. 
Moreover, the construction developed in j4| for the said purpose makes use of 
the curves Vq^ corresponding to the pressures don, which are wrongly regarded 
as curves with isolated points of contact as it was noted and discussed above. 

Below, an alternative approach is presented for constructing equilibrium 
ring (tube) shapes with lines (areas) of contact based on the same "similarity" 
idea that actually arises out of the following property of Eqs. ( 122|) and (!23|) . 

Under the transformation (s, k) i — )■ (s/A, Ak), where A is an arbitrary real 
number, each equation of form (1221) corresponding to certain constants /i and 
(J transforms into an equation of the same form but with new coefficients: 
/i I — )■ X^fi, a I — )■ The same holds true for equation f l23|) if e i — )■ in 
addition. In other words, equations ( 122|) and fl23|) are invariant with respect 
to the similarity transformation A : (s, k\ fi, a, e) i — > (s/A, Ak; A^yU, A^cr, X^e). 
Consequently, the parametric equations (125|) imply that the shapes whose 
parameters are related by such a transformation A are similar, the respective 
scaling factor being 1/A. Accordingly, if a closed curve F is scaled in this 
way, then its length L and area A change to L/X and A/X"^, respectively. 
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Thus, given n > 2, let the curve Ten of length Lc„ = 27r be the equilibrium 
shape with points of contact corresponding to the contact pressure and 
let r be the shape (of the same length) with lines of contact corresponding 
to a pressure a > acn- The curve f is constructed in two steps. First, scaling 
the curve Fen with a factor {fTlacn)^^^ one obtains another curve Fen which 
has the same number of contact points because it is similar to the curve Fen 
but corresponds to the pressure a and its length is Lcn — '2-Ti{(Tcn/ (^Y^^ < Lcn- 
Then, the curve f is obtained by substituting each point of contact of the 
curve Fen by a line segment of length 27r(l — ((7cn/'5')^'^^)/n along the respective 
symmetry axis of the curve Fen so as its total length to become 27r. 

Examples of shapes with hues of contact are presented in Figs. 7 and 8. 




Figure 7: Shapes with lines of contact corresponding to: (a) a = 10.34 (2-fold symmetry); 
(b) a = 81.81 (3-fold symmetry); (c) a = 207.2 (4-fold symmetry). 




Figure 8: Shapes with lines of contact corresponding to: (a) a = 400 (6-fold symmetry); 
(b) cr = 800 (9-fold symmetry); (c) a = 1500 (12-fold symmetry). 
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It is clear that the tangent, normal and position vectors of a shape f with 
lines of contact constructed in the foregoing way are continuous at each point 
of the curve F. However, its curvature suffers jumps at the end points of the 
line segments used to substitute the contact points of the respective auxiliary 
curve Ten because the limit values of the curvature from the bent parts of 
the curve and from the line segments are — -\/2/I 7^ and zero, respectively. 
Consequently, the moment and force also suffer jumps at the aforementioned 
points since their limit values from the bent parts of the curve are 



Mfe = (v^+/€°) , Nb = 0, Qfe = ±Z^^p(-v^), (65) 

while along each line of contact the resultant pressure is zero and 

Mi = -Dk°, Ni = 0, Qi = 0. (66) 

Equations (!65|) and (!66|) follow by the constitutive equation (120|) . the general 
solution ([2ID of Eqs. ([7D - © and Eqn. (|23l). 

Thus, the local balances and (E]) of the force and moment are violated 
for the shapes with lines of contact. Fortunately, however, the total balances 

^F'{s)ds = -^pn{s)ds, (67) 

^M'{s)ds = -^F{s)-n{s)ds, (68) 

of these quantities are satisfied. Indeed, Eqs. (E7|) and (1^ hold on the curve 
Ten since it corresponds to an equihbrium shape without jump discontinuities 
of the force and moment. On the other hand, F = Fen U {Line segments} 
and the integrals in Eqs. fl67|) and fl68|) taken along the line segments are 
equal to zero because here p = 0, M'(s) = and F(s) = 0. 
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In our opinion, this property of the constructed curves f allows this shapes 
to be regarded as equilibrium ring (tube) shapes with lines (areas) of contact 
at least in the week sense discussed above. 

In the light of the results presented in this Section, it should be remarked 
that the similarity law (5.4) obtained in ^ Section 5], which concerns the 
conductivity of a buckled tube conveying an incompressible viscous fluid, has 
to be revised. Actually, this law, which expresses the conductivity of a tube 
with areas of contact through that of a tube whose cross sections have just 
points of contact, should be replaced by the following one 



(T,„ \ 4/3 



a 



C{a) = {^] C(ae„), a>ae„, (69) 



where C{a) denotes the conductivity of a tube subject to pressure a. It is 
necessary to do so because (Ten is the unique pressure for which there exists 
an equilibrium tube shape of n-fold symmetry whose cross sections have only 
isolated points of contact. 

9. Concluding remarks 

In the present paper, the problem for determination of the equilibrium 
shapes of a circular inextensible elastic ring (tube) subject to a uniform 
hydrostatic pressure is reexamined. For the first time, more than a centiiry 
ago, this problem was stated and studied by Maurice Levy in his memoir 7|. 

Here, a concise derivation of the most important facts established in 7|, 
see 1° - 3° (p. 9), concerning the existence of a "centre of the elastic forces", 
following by Eq. ffTOj) . and the properties reflected by Eqs. f lTTj) and f lT2|) is 
given in Section 3. Then, the parametric equations of the equilibrium shapes 
are expressed through the forces and slope angle, see Eqs. ffT7|) . 
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It is noteworthy that neither the relations fllUp - (fT^ nor the forms of 
the parametric equations (|T71) or boundary conditions f|T8|) and f|T9|) depend 
on the particular material properties or stress-free configuration of the ring. 

A concise justification of the symmetry of the ring shapes discussed by 
many authors (usually without going into much detail) is given in Section 4. 
It is shown analytically that each ring shape is symmetric and can be obtained 
by successive reflections of its part corresponding to the first half period of 
the curvature with respect to its symmetry axes. 

Further, assuming that the stress-free configuration of the ring is a circle 
of radius p, the case in which the linear constitutive equation f l20|) holds is 
considered. In this case, the equilibrium state of the ring is determined by the 
periodic solutions of the nonlinear ordinary differential equation (^E^ for the 
ring curvature, which are such that the closure condition ( 127|) holds. In fact, 
the shape of the ring is determined explicitly by the parametric equations 
f l25p and the values of the moment and forces acting along the ring axes are 
given by Eqs. fl20|) and f l2T|) . respectively. Explicit analytic expressions for 
all periodic solutions of Eq. fl23|) and for the corresponding slope angles are 
presented in Section 5. 

In Section 6, the determination of the equilibrium shapes corresponding 
to a given pressure a is reduced to the computation of the common solutions 
of two transcendental equations ( H2|) and (143|) . 

Shapes with isolated points of contact are studied in Section 7. It is 
shown that the pressures at which such a shape is attained can be obtained 
computing the common solutions of Eqs. fj42|) . f H3|) and flM|) in the case 
n = 2, or Eqs. (02]), (03]) and ([60]) for n > 2. 
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The most important results achieved in Sections 6 and 7 follows. 
In contrast to the assertion in ^ that for each mode n there exists a range 
of pressures for which the respective ring shape has only isolated points of 
contact, we found, solving numerically Eqs. fH2|) . fH3|) and f lM|l or fl60l) . that 
for each mode 2 < n < 15 there is a unique such pressure, namely acn- 
Moreover, for all modes in the range 2 < n < 15 our computations based on 
the procedure described in Section 6 show that if the applied pressure a is 
such that < cr < Ccn, then the corresponding buckled shape of n mode 
is simple, while for a > acn this shape always has points of self-intersection. 
Our conjecture is that this behaviour is inherent to all modes. 

Section 8 concerns the equilibrium ring (tube) shapes with lines (areas) of 
contact that are expected to occur for pressures greater than the respective 
contact pressure instead of the self-intersecting shapes (unnatural for tubes 
and planar rings) predicted by the considered model. Here, the construction 
of these shapes is based on the similarity properties of Eqs. and 
following in general outline the idea suggested by Flaherty et al. in [4^, but 
the uniqueness of the contact pressures is taken into account. The shapes 
obtained in this way are shown to satisfy the total balances f l67|) and f l68|) of 
the respective forces and moments, which is good reason to consider them 
as equilibrium shapes. Finally, the expression for the so-called similarity law 
obtained in j4| is revised, see Eq. (169|1 . 

The interested reader can find the Mathematica® notebooks developed 
for the solution of the transcendental equations of Sections 6 and 7 as well 
as the notebooks developed for the construction of the shapes described in 
Section 8 at http://www.bio21.bas.bg/ibf/dpb_files/mfiles/. 
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